d-Wacve Pairing in the Presence of Long-Range Coulomb Interactions 
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The one-band extended Hubbard model in two dimensions near band-filling i is solved in the 
fluctuation exchange approximation, including the long-range (1/r) part of the Coulomb interaction, 
up to 4th neighbor distance. Our results suggest that d^2_y2 pairing in the Hubbard model is 
robust against the inclusion of long-range Coulomb interactions with moderate 1st neighbor repulsion 
strength Vi. d^2_y2 pairing is suppressed only at large Vi (> 0.25(7-0.4(7), due to incipient charge 
density wave instabilities. 
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Several cuprate superconductors appear to exhibit a 
pairing state of d^2_y2 symmetry [Q. Such a non-s- 
wave state was also found early on in microscopic pairing 
theories, based on antiferromagnctic (AF) spin correla- 
tions in the two-dimensional (2D) Hubbard model [^. 
While incorporating the local on-site U Coulomb repul- 
sion, the Hubbard-based d^2_y2 pairing theories have so 
far largely ignored the effects of the spatially more ex- 
tended Coulomb matrix elements. Although unaffected 
by the on-site U, non-s-wave pairing states can, in princi- 
ple, be severely suppressed by the longer range Coulomb 
repulsion, due to the spatially extended nature of their 
Cooper pair wavefunctions. For any pairing theory it is 
therefore of central importance to establish whether non- 
s-wave pairing survives the extended repulsion. 

Here, we address this question in a self-consistent dia- 
grammatic framework, the fluctuation exchange (FLEX) 
approximation applied to an extended Hubbard 
model. We show that dx2-y2 pairing is (i) robust 
against extended Coulomb terms of realistic strength and 
(ii) suppressed only if the extended repulsion becomes so 
strong that it induces charge density wave (CDW), rather 
than AF spin density wave (SDW) instabilities in the ^- 
filled Hubbard system. This CDW scenario is apparently 
not realized in the undoped cuprates . 

We start from the 2D extended Hubbard Hamiltonian, 
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U/t^ 8-12 Absorbing dielectric screening effects 

from the insulating non-Hubbard electron background 
into an effective Vi , one estimates Vi/t^ 2-3 [|| . Addi- 
tional screening from phonons reduces this to Vi/t<Q.3- 
0.5 but should really be treated expHcitly in a more 
realistic electron-phonon model, because of the low (i.e., 
potentially relevant) phonon energy scale. 

In Fig. |l|(a), we show typical (3rd order) FLEX con- 
tributions to the single-particle self-energy S, written in 
terms of the bare particle-hole vertices Vd for density and 
Vm for magnetic fluctuation exchange. In Fig. |l|(b) and 
(c) , we show the bare vertex Vg and, respectively, the 1st 
and typical higher (3rd) order contributions to the renor- 
malized vertex Fg for the singlet particle-particle interac- 
tion. Wc do not include particle-particle fluctuation ex- 
change in E. Hence, within our FLEX approximation 
no Aslamazov-Larkin type diagrams contribute to Fg 
at or above Tc (see also previous applications in Ref. |0 ) . 

The finite cutoff, V^(AR) = for |AR| > r„, makes 
a numerical solution of the FLEX equations |^ feasible. 
The key algorithmic ingredients are (i) introduction of a 
mixed real-space and momentum-space basis set [|| for 
two-body propagators; and (ii) application of a numerical 
renormalization group to calculate Matsubara frequency 
sums (l^]. The vertices (Vd, Vm, Vg) are MxM matrices, 
where M is the number of AR with y(AR) ^ 0. For 
example, in the 1st neighbor (n = 1) Coulomb model, 
M — 5 and their non-zero matrix elements ai Q = (Q, il) 
[see Fig. |^(a) and (b)] are 



where creates a hole with spin a at site Ri on an 
N = Lx L square lattice with periodic boundary condi- 
tions, and riia- — c\^Cia. In tij, we include only a 1st 
neighbor hybridization t and the chemical potential /i. 
The Coulomb matrix element Vij = V(R.i — R^) com- 
prises the Hubbard on-site repulsion U = V{0) and an ex- 
tended l/r-part, V{AR) = Vi/|AR| for < |AR| < r„, 
up to the 7ith neighbor distance r„. Here Vi denotes 
the strength of the 1st neighbor repulsion and AR is in 
units of the lattice constant a with a = 1 in the follow- 
ing. For the cuprates, we assume t ~ 0.3-0.5 eV and 



Fd(Q;AR; AR') - 

C/ + 4Fi(cosQ^ +cos(5j^), AR = AR' = 
-Vi, 

K.(Q;AR; AR') = 

-u , 

-Vi , AR = AR' = ±x, ±y 

Fs(Q;AR; AR') = 

U , AR = AR' = 

Vi/2, AR=±AR' = ±x,±y. 



AR = AR' ±x,±y 
AR = AR' = 
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FIG. 1. Diagrams included in the FLEX approxi- 
mation: (a) Representative higher (3rd) order contribu- 
tion to the fluctuation self-energy, E(fc). (b) Bare sin- 
glet pairing vertex, Ki(Qs; AR; AR') with total pair mo- 
mentum-energy Qs = (Qs,r2s). (c) 1st and typical higher 
(3rd) order contributions to the full singlet pairing vertex, 
rs(Qs; AR,tj;AR',t^')- 



The superconducting transition occurs when = 1 
where is the largest eigenvalue of the integral kernel 
for the particle-particle ladder series with Fg as the inter- 
action vertex |^ at total pair momentum-energy (5s = 0. 
The Tc results reported here are for 16 x 16 lattices and 
agree with 32x32 lattice results to within 5% at selected 
parameter values we have checked. 

In Fig. 1^ we plot the superconducting transition tem- 
perature Tc as a function of Vi for different hole con- 
centrations and different U in the 1st neighbor [n = 1) 
and 2nd neighbor [n = 2) Coulomb repulsion models. 
As in previous FLEX calculations for Hubbard-like mod- 
els, the symmetry of the dominant pairing instability is 
dx^-y^ P,p|, pO| , pT| . The striking feature in Fig. |^ is the 
remarkable insensitivity of to Vi over a wide param- 
eter range. An eventual abrupt down-turn in Tc occurs 
only when Vi approaches the crossover l/™*^ from SDW 
to CDW fluctuation behavior, discussed below. 

To understand this behavior, we analyze the decompo- 
sition of Fg in Fig. |l|(c) for the case of the 1st neighbor 
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FIG. 2. d^2_y2 transition temperature Tc vs. 1st neighbor 
Coulomb repulsion, Vi, at (a) U = 8t and (b) U = At for 
hole densities (n) = 1.12 (circles), 1.16 (squares), 1.20 (dia- 
monds), and 1.24 (triangles). Open symbols are for the 1st 
neighbor (n = 1), closed symbols for the 2nd neighbor (n = 2) 
Coulomb model. Vertical dashed (dot-dashed) lines mark the 
crossover V{^'^ from SDW to CDW fluctuation behavior in the 
1st neighbor (2nd neighbor) model, as described in text. 



(ri = l) model. We consider T ^{0; AR^inT; AR, inT) as 
a measure of the effective interaction between holes sep- 
arated by AR in the Qg — singlet pairing channel at 
temperature T. Note that Vi affects Tc via the 1st order 
(Hartree-Fock) contribution to E (not shown in Fig. |^), 
via the 1st order contribution [Fig- 0(c)] to Fg, and 
via the charge and spin fluctuation contributions, that 
is, via the Vd- and V^n-ladder series which enter into E 
[Fig. 1(a)] and Fg [Fig. g(c)]. 

The primary reason for the robustness of the d-wave 
pairing against moderate 1st neighbor repulsions, Vi ^ 
Vf^'^, is simply the raw strength of the effective electron- 
electron interaction generated by the spin fluctuations, 
primarily with momentum transfer Q [see Fig. |l|(a)] near 
Q = (7r,7r). For example, at 1st neighbor distance, U = 8t 
and hole density (n) = (^^rij^) = 1.16, this attractive 
contribution to Fg is of order 32t. This clearly dominates 
over extended Coulomb effects at moderate Vi. 

At U = 8t, the primary effect of a moderate Vi on Fg is 
a slight enhancement of the 1st neighbor attraction, by 
net amounts ^ 0.2t. This additional VL-induced attrac- 
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tion is caused by Q ~ charge and by Q '--^ (tt, tt) spin 
fluctuation contributions to Fg, which overscreen the di- 
rect 1st order repulsive Fg contribution, Vs, in Fig. ||(c). 
The reason for the sUght suppression of Tc by a moder- 
ate Vi at U = 8t is therefore not the Vi -effect on Fs, but 
rather its effect on E. By increasing S, Vi suppresses 
the single-particle spectral weight near the Fermi energy 
which, in turn, tends to reduce Tc. This S-effect out- 
weighs the Fg-effect of Vi . 

As indicated by the U — -it results in Fig. ^(b), a large 
f/- value is important for stabilizing the (1^2 _y2 pairing 
against extended Coulomb effects. To understand this, 
note first that the spin fluctuation contributions to Fg 
and E are strongly [/-dependent. For example, at U = 
4t, Vi = and {n) = 1.16, the 1st neighbor attraction 
induced by the spin fluctuations in Fg is reduced to 3.47t, 
compared to ~ 32t in the U = 8t case. Despite this 
order of magnitude decrease, Tc for U — At is slightly 
higher, due to the accompanying reduction of the self- 
energy effects discussed above. However, for this U = 4:t 
parameter set, Tc is also suppressed much more rapidly 
by turning on Vi. This is due to the fact that the spin 
fluctuation contribution to the 1st neighbor attraction in 
Fg is now suppressed, rather than enhanced, by Vi. As 
a consequence, the net effect of Vi on Fg is to suppress 
the 1st neighbor attraction, for example by about 0.6t for 
Vi = 0.5Vi^'^. Due to the reduced overall strength of Fg, 
this is a relatively much larger effect than in the U — 8t 
case and therefore has a much larger effect on Tc. 

The SDW-CDW crossover V{^'^ in Fig. | is opera- 
tionally defined as that Vi -value where the maximal 
"density" eigenvalue A'^ becomes equal to the maximal 
"magnetic" eigenvalue A™. Here, A^ and A™ denote the 
largest eigenvalues of the integral kernels for the Vd- and 
VJn-based particle-hole ladder series entering into E and 
Fg, shown in Fig. |](a) and (c). This criterion for V™'^ 
may give only rough estimates of the actual physical 
crossover since, instead of the full interaction vertices Fd 
and Fni it is obtained from the bare interaction ver- 
tices Vd and Vra- The estimates are reasonable vis-a-vis 
t^Q limit results discussed below. 

As Vi approaches VY"*^, CDW fiuctuations rapidly take 
over and due to their nearly singular nature preclude 
convergence of our self-energy iteration for T ^ Q.\t 
and Vi > F™^. In the 1st neighbor {n = 1) Coulomb 
model, the CDW fluctuations exhibit maximum A*^ near 
Q = (tt, tt). They produce a repulsive 1st neighbor contri- 
bution to Fg and, via E, they weaken the strength of the 
the attractive Q~ (tt, tt) spin fluctuation contribution to 
Fg when Vi V{^'^. The Q ~ (tt, tt) charge fluctuations 
also generate an on-site attraction which is too weak, 
however, to overcome the on-site repulsion due to U and 
Q ~ (tt, tt) spin fluctuations. No pairing instabilities are 
found down to T^^ O.lt for Vi > V™'^; the dominant at- 
tractive symmetries are odd-w s-wave in the triplet and 
even-w s-wave in the singlet channel with respective pair- 



ing eigenvalues A' 0.35-0.45 and A" ~ 0.25-0.4 in the 
10-20% doping range. d^2_y2 pairing is suppressed. The 
singlet s-wave attraction suggests that CDW fluctuations 
may enhance conventional phonon-mediated pairing. In 
the 2nd neighbor (n = 2) Coulomb model, the SDW- 
CDW crossover is similar, but driven by CDW fiuctua- 
tions with maximum A'* near Q = (0, tt) and (7r,0). 

The physical origin of the SDW-CDW crossover can 
be understood by considering the exactly solvable ionic 
limit, i— >0, at ^-filling ||^. Large U favors a groundstate 
with single occupancy at each site and AF SDW order for 
\t\ =0+. Large Vi, in the 1st neighbor Coulomb model, 
favors a CDW groundstate with modulation wavevector 
Q* = (tt, tt) and alternating double and zero occupancy in 
the respective two sub-lattices. The transition between 
these two states occurs when Vi reaches ||^ y-md|t-o ^ 
\U which coincides with the exact crossover V™*^ — \U 
shown in Fig. || In the 2nd neighbor Coulomb model, 
the competing CDW state has Q* = (0, tt) or (tt, 0), with 
alternating doubly and zero occupied rows or columns, 
and the transition occurs at V™'^'*~° = -\/2[//4 = 0.354C/. 
This is roughly consistent with the approximately T- and 
(ri)-independent crossover 

ymdg^Q 406[/ from Fig. |(a). 
This analysis suggests that the crossover is driven by in- 
cipient CDW instabilities in the i-fiUed Mott-Hubbard 
insulator. We emphasize that such CDW instabilities 
and their effects on pairing cannot be treated in strong- 
coupling versions of the Hubbard model, such as the t-J 
model CDW fiuctuations in the ^-filled system are 
precluded from the outset by the no-double-occupancy 
Hilbert space constraints imposed in such models. 

By exploiting the retarded nature of the pairing po- 
tential, the pair wavefunction can evade the destruc- 
tive effects of the instantaneous Coulomb repulsion , 
provided the characteristic frequency Ob of the pairing- 
mediating boson is much lower than characteristic elec- 
tronic energy scales such as the Fermi energy ep. This 
Coulomb "pseudopotential" effect |l^ is crucial for 
conventional phonon-mediated superconductors, where 
ile/ep ^ 10^^-10^"^. However, in the present problem, 
ile/eF ^ 10~^ is not very small, since the spin fiuctua- 
tion frequencies Ob extend up to a sizeable fraction of the 
electronic bandwidth ||l^,|l^,^. By a variational analy- 
sis, based on restricted trial pair wavefunctions, we can 
indeed show that the pseudopotential effect raises \^ and 
Tc by no more than a few percent, i.e., it is largely inop- 
erative in spin fiuctuation mediated pairing. 

Next, we consider longer-range 1/r Coulomb terms. As 
shown in Fig. ^(a), the inclusion of 2nd neighbor (n = 2) 
Coulomb terms suppresses Tc more strongly at low Vi 
than in the 1st neighbor model. However, the d-wave 
pairing also survives to larger Vi, since Vj™'^ is increased 
in the 2nd neighbor model. Comparable Vj™'^ are found 
in longer-range models with n>2. 

As shown in Fig. |[ the dependence of Tc on the cutoff 
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der in undoped insulating cuprates appears to rule out 
the latter scenario. 
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FIG. 3. d^2_y2 Tc vs. cutoff radius r„ of the extended 
Coulomb interaction. Here, r„ = 0, 1, \/2, 2, and for 
n = 0, 1, 2, 3, and 4th neighbor cutoff, respectively. Open 
symbols are for U/t — 8, Vi/t — I; closed symbols are for 
U/t = 4, Vi/t = 0.5; with (n) = 1.12 (circles), 1.16 (squares), 
1.20 (diamonds), and 1.24 (triangles). 



radius r„ at fixed Vi rapidly saturates for n>2. This can 
be understood by noting that the long-range {n>3) tail 
of the extended Coulomb potential is strongly suppressed 
by the metallic screening, while the the d-wave pairing 
wavefunction is mostly concentrated at 1st neighbor dis- 
tance. Hence the effect of the additional {n > 2) direct 
repulsion terms on the pairing is rapidly reduced with in- 
creasing r„. We are thus confident that our conclusions 
based on finite-range (n < 4) models are not fundamen- 
tally altered in the infinite-range 1/r model. 

For comparison, we have also explored the effects of an 
added extended Coulomb repulsion on the c?-wave Tc in 
recently proposed phenomenological spin fluctuation ex- 
change models 



The d-wave pairing is much more 
rapidly suppressed with Vi in such models, primarily due 
to the lack of screening, i.e., due to the fact that the spin 
fluctuation mediated pairing potential, as extracted from 
fits to experimental data, is "rigid" and does not become 
modified by the extended Coulomb repulsion. 

In summary, we find that, at the level of the FLEX ap- 
proximation to the one-band Hubbard model, spin fluc- 
tuation mediated (1^2 _y2 pairing is robust against ex- 
tended 1/r electron-electron Coulomb repulsions of re- 
alistic strengths. The extended part of the Coulomb 
interaction does not cause a significant change in the 
spin fiuctuation spectrum nor in the pairing vertex of 
the Hubbard model for moderate strengths Vi. The ro- 
bustness of the model against the inclusion of extended 
Coulomb terms increases with increasing on-site repul- 
sion U. With increasing real-space cutoff radius r„, the 
extended Coulomb effects are almost entirely saturated 
after the inclusion of the 2nd neighbor repulsion terms. A 
strong suppression of (1^2 _y2 pairing is found only when 
the extended Coulomb repulsion strength Vi becomes so 
large that it induces CDW instabilities in the ^-filled 
Mott-Hubbard insulator. The ubiquity of AF SDW or- 
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